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Abstract —We provide nonasymptotic upper and lower bounds 
on the sum-rate capacity of Rayleigh block-fading multiple-access 
channels for the set up where a priori channel state information is 
not available. The upper bound relies on a dual formula for chan¬ 
nel capacity and on the assumption that the users can cooperate 
perfectly. The lower bound is derived assuming a noncooperative 
scenario where each user employs unitary space-time modulation 
(independently from the other users). Numerical results show that 
the gap between the upper and the lower bound is small already 
at moderate SNR values. This suggests that the sum-rate capac¬ 
ity gains obtainable through user cooperation are minimal for the 
scenarios considered in the paper. 


I. Introduction 

The multiple-access channel (MAC) models a scenario where 
two or more noncooperating users communicate with a single 
receiver. This scenario is relevant for the uplink of wireless 
cellular networks, where the users may be mobile terminals and 
the receiver may be a cellular base station. In this paper, we 
consider the setup where neither the users nor the receiver have 
a priori information on the realization of the fading process. 
Such a situation arises in high mobility scenarios where it is not 
desirable for the receiver to feed back channel state information 
(CSI) to the users because it may be outdated [1]. It may also 
arise in the initialization phase of a communication link, e.g., 
when a mobile terminal joins a cellular network. Capacity anal¬ 
yses under assumption of no a priori CSI have the advantage of 
capturing the cost of estimating the fading channel, and, hence, 
yield more realistic throughput estimates than analyses based 
on the assumptions of perfect CSI [2]—[4]. Indeed, under the 
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assumption of no feedback from the receiver to the user termi¬ 
nals, these capacity analyses provide a fundamental limit on the 
performance of every communication scheme, irrespectively of 
whether it relies on explicit channel estimation or not. 

In this paper, we consider the case where the users as well as 
the receiver are equipped with one or more antennas. We shall 
focus on the so-called Rayleigh block-fading model [5], [6]. The 
two key features of this model are that i) the fading coefficients 
associated to the channels between each transmit and receive 
antenna pair are independent and identically distributed (i.i.d.) 
circularly symmetric complex Gaussian random variables; ii) 
each fading coefficient remains constant over t c channel uses 
before changing to a new independent realization. The parame¬ 
ter t c , which is the ratio between the channel coherence time 
and the symbol duration, will be referred to in this paper as 
coherence interval. 

The capacity of Rayleigh block-fading channels under the 
assumption of no a priori CSI has been studied extensively in 
the point-to-point case. Specifically, considering a system with 
n t transmit and n r receive antennas, Marzetta and Hochwald 
proved that the capacity-achieving input matrix X £ C n,xtc 
can be expressed as [6, Thm. 2] 

X = DQ (1) 

where D is a real n t x n t diagonal matrix whose diagonal entries 
have a joint probability density function (p.d.f.) that is invariant 
to permutation of its arguments, and Q £ C n,xtc . independent 
of D, is an isotropically distributed matrix with orthonormal 
rows (truncated unitary matrix). 

In spite of the partial characterization of the capacity- 
achieving input distribution provided in [6, Thm. 2], no closed 
form expression for capacity is available to date. However, 
the high-SNR capacity behavior is well understood. Indeed, 
extending a result obtained for the single-input single-output 
case in [5], Zheng and Tse [7] proved that in the high-SNR 
regime, the capacity C of a n t x n r multiple-input multiple- 
output (MIMO) Rayleigh block-fading channel behaves as 

C(p) = n* (1 - log(p) + 0{ 1). ( 2 ) 

Here, p stands for the SNR, n* = min(n t , n r , [_i c /2j), and 0(1) 
indicates a function whose magnitude is upper-bounded by a 
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finite constant for sufficiently large SNR values. The asymptotic 
expression (2) can be tightened to [7], [8] 

C{p) = n* ^1 - ^ log(p) + c + o(l) (3) 

where c is the constant given in [8, Eq. (9)] and o(l) denotes 
a function that vanishes as p — > oo. For the case when t c > 
n t + n r . one can achieve (3) by choosing D in (1) to be a 
scaled identity matrix [7]. The resulting probability distribution 
on X is commonly referred to as unitary space-time modulation 
(USTM) [5]. Whenf c < n t +n r , the matrix D must be chosen so 
that its diagonal entries are distributed as the square root of the 
eigenvalues of a Beta-distributed random matrix of appropriate 
size [8]. The resulting probability distribution on X is referred 
to in [8] as Beta-variate space-time modulation (BSTM). 

Nonasymptotic (i.e., finite-SNR) lower bounds on the capac¬ 
ity of point-to-point MIMO Rayleigh block-fading channels 
have been obtained for specific probability distributions on X. 
Specifically, an i.i.d. Gaussian input distribution is considered 
in [9], USTM in [10], and BSTM in [11], where the analysis is 
also extended to Rician block-fading and to land mobile satellite 
channels. A key tool in the derivation of these nonasymptotic 
lower bounds is the Itzykson-Zuber integral [12, Eq. (3.4)], 
which allows one to obtain a closed-form expression for the 
conditional probability distribution of the channel output given 
the diagonal input matrix D in (1). The method employed so far 
to assess the tightness of the bounds obtained in [9]—[11] is to 
compare these lower bounds with the asymptotic expansion (3) 
(with the o(l) term omitted). Unfortunately, this method is not 
conclusive because the error incurred by omitting the o(l) term 
in (3) is not quantified. A simple capacity upper bound can be 
obtained by assuming that a genie provides the receiver with 
perfect CSI. However, this bound is tight only when t c is large 
and the channel estimation overhead negligible. 

Leaving the point-to-point case and moving to the MAC, we 
note that the independence constraint on the signals transmitted 
by the various users implies that the partial characterization of 
the capacity-achieving input distribution obtained for the point- 
to-point case in [6, Thm. 2], as well as the asymptotic capacity 
expansions in [7], [8] and the nonasymptotic capacity lower 
bounds in [10], [11] do not carry over to the MAC sum-rate 
capacity. The only exception is the i.i.d. Gaussian lower bound 
obtained in [9], which also applies to the MAC because the trans¬ 
mission of i.i.d. Gaussian signals does not require coordination 
among the users. 

Coarse upper and lower bounds on the MAC sum-rate ca¬ 
pacity for the case when n u single-antenna users communicate 
with a receiver equipped with n r antennas are provided in [13]. 
By examining these bounds in two different asymptotic regimes 
(high SNR and large t c . for a fixed n u /t c ratio) the authors con¬ 
jecture that the sum-rate capacity is maximized when n u = t c . 
In the same paper, the authors pose the question of whether 
the constraint that the users transmit independent signals yields 
a sum-rate capacity prelog 1 that is strictly lower than the one 
achievable when the users can cooperate perfectly, and the MAC 

1 The prelog (a.k.a. multiplexing gain) is the asymptotic ratio between the 
(sum-rate) capacity and log(p) in the limit p —>■ oo. 



Fig. 1. Upper and lower bounds on the MAC sum-rate capacity: 4 single¬ 
antenna users; receiver with 4 antennas; coherence interval of 10 channel uses. 

reduces to a point-to-point MIMO channel, for which (2) holds. 
It follows from [7, Sec. V] that this is not the case, provided 
that the users are able to transmit orthogonal pilot signals, used 
at the receiver to estimate the channel. 

Lin and Moser [14] characterized the high-SNR behavior of 
the sum-rate capacity of an i.i.d. Rician-fading MAC (block¬ 
fading channel with coherence interval t c = 1). They showed 
that the sum-rate capacity grows double-logarithmically in SNR, 
and that the sum-rate capacity maximizing strategy at high SNR 
is to switch off all users but one. 

Contributions: We present nonasymptotic (i.e., finite-SNR) 
upper and lower bounds on the sum-rate capacity of Rayleigh 
block-fading MACs. Similarly to the nonasymptotic capacity 
lower bounds previously reported in [9]—[11], our bounds are 
not in closed form, but they can be evaluated numerically. Our 
upper bound is obtained by assuming that the users can perfectly 
cooperate, which turns the MAC into an equivalent MIMO point- 
to-point channel. In addition, we use the duality upper bound on 
mutual information reported in [15, Eq. (186)]. As auxiliary out¬ 
put distribution in the duality step, we choose the one induced 
by USTM inputs in the absence of additive noise. This method 
was used in [8] to establish the asymptotic expansion (3) for 
the case n t < min(n r , [f c /2j). Here, we provide a finite-SNR 
analysis, which generalizes to MIMO the one reported in [4] for 
the single-input single-output (SISO) case. 

The nonasymptotic lower bound is obtained by allowing each 
user to transmit a USTM signal over the available antennas. 
We call the resulting input distribution MAC-USTM. Note that 
MAC-USTM does not yield a global USTM input distribution, 
because independence among users imply that orthogonality 
among them cannot be enforced. Numerical evidence (see Fig. 1 
and Fig. 2) suggests that for the scenarios considered in the 
paper: 

• Our upper and lower bounds characterize accurately the 
sum-rate capacity (see Fig. 1). 

• The gain in sum-rate capacity obtainable by allowing user 
co-operation is minimal. This follows because the upper 
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Fig. 2. Upper bound, MAC-USTM lower bound, and i.i.d. Gaussian lower 
bound on the MAC sum-rate capacity: 4 single-antenna users; receiver with 4 
antennas; coherence interval of 10 channel uses. 


bound is obtained under the assumption of perfect coopera¬ 
tion between users. A similar observation has been recently 
reported in [16] for the nonfading asynchronous Gaussian 
MAC. 

• Orthogonal pilot transmission is not required to obtain 
rates close to the full-cooperation case (cf. [13], [7, Sec. 
V]). 

• The MAC-USTM lower bound is tighter than the i.i.d. 
Gaussian lower bound, although the gain is marginal when 
the users are equipped with a single antenna (see Fig. 2). 
This suggests that i.i.d. Gaussian inputs are almost sum- 
rate capacity optimal already at moderate SNR values, 
which confirms an observation reported in [9] for the point- 
to-point MIMO case. 

Notation: Uppercase letters denote matrices, 2 lowercase let¬ 
ters designate scalars, and boldface letters denote random quan¬ 
tities. The superscript f stands for Hermitian transposition. For 
a random variable x with p.d.f. f x (x), we write x ~ / X (x). 
When two random variables x and y have the same p.d.f., we 
write x ~ y. For a matrix X , we let X, :j denote its entry on the 
?th row and the jth column. With detwe indicate the 
determinant of a matrix whose entry on the 2 'th row and the jth 
column is given by for some arbitrary function /(•, •). 

It will turn out convenient to define the following function 


T(n) 


n r «> 

i— 1 

i, 


n £ N 


n = 0 


(4) 


where T(-) is the Gamma function [17, Sec. 6.1]. With /?(-, •) 
we denote the Beta function [17, Sec. 6.2] and ip(-) denotes 
the Digamma function [17, Sec. 6.3]. The identity matrix of 
dimension nxn is denoted by I n . We let diagjdi, d, 2 ,... ,d n } 


2 We do not distinguish between vectors and matrices. We treat n- 
dimensional vectors as 1 x n or n x 1 matrices. 


be the diagonal matrix with entries d\ , d ,2 , ■ ■ •, d n on its diago¬ 
nal. For a Hermitian matrix A £ C" x ” with ordered eigenval¬ 
ues ai > a 2 > ■ ■ ■ > a n we denote the determinant of the 
Vandermonde matrix constructed from a±, < 22 ,..., a n as 

V(A)= n (a* — a j). (5) 

We shall often use the following two functions 

/C(A, k) = V(A) det(A) fe (6) 

where A £ C nxn is Hermitian and k £ N, and 

n—1 k 

7 (x,n) =e x ~Y,^ (7) 

fc=0 

where i£l and n £ N. The set of unitary matrices in <C nxn 
is denoted by U(n) ( unitary group) and the set of matrices 
U £ C" xm ,m > n with UW = I n is denoted by S(n,m ) 
(Stiefel manifold). With E x [/(x)], we denote the expectation of 
the function /(x) over the random variable x. We let h(x.) de¬ 
note the differential entropy of a continuous random variable x; 
furthermore, X(y; z) stands for the mutual information between 
the random variables y and z. The set of all diagonal matrices in 
JR" x m with ordered and distinct positive entries on their main 
diagonal is denoted by U" x m . With 0„ x m we indicate the nxm 
zero matrix. We use CA/”(0, o 2 ) to denote a circularly symmetric 
complex Gaussian random variable with zero mean and variance 
cr 2 , and Beta(a, b) to denote a Beta-distributed random variable 
with parameters a and b. 

II. System Model 

We consider a Rayleigh block-fading MAC where n u users 
communicate with a receiver having n T antennas, and the chan¬ 
nel coherence interval is t c (same for all users, which corre¬ 
sponds to a scenario where users with similar mobility require¬ 
ments are scheduled together). We assume that each user is 
equipped with one or more antennas and denote by n,; the 
number of antennas at user i. i = 1,..., n„. The received 
signal Y £ C n,xtc over an arbitrary coherence interval can be 
compactly written in matrix notation as follows: 

n u 

Y = ^S l X,+W. (8) 

»=1 

Here, X,; £ C n ' xtc denotes the signal transmitted by user i over 
the coherence interval, and the matrix S, £ C n ' x r '"‘ contains 
the fading coefficients associated to the channels between each 
transmit antenna of user i and the receive antennas, within the 
coherence interval. We assume that Sj has i.i.d. CA'jO. 1) entries 
and that the channel matrices {Sj}"2. 1 are independent. Finally, 
the matrix W £ C" rXtc , whose entries are i.i.d. CAf(0, 1)- 
distributed, denotes the additive noise. Let 

n u 

n t = ^2rii (9) 

*= 1 

be the total number of transmit antennas. We can rewrite (8) as 
Y = SX + W (10) 
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where 


Here, 


S n „] er' xn ' (li) 


G C n,xtc . (12) 

We assume that W and S are independent, and that their proba¬ 
bility law does not depend on X. Throughout the paper, we also 
assume that t c > max(n t , n r ) and focus on the scenario where 
neither the transmitter nor the receiver have prior knowledge of 
matrix S (no a priori CSI). It will turn out convenient to define 
the following two constants 

( = min(n t ,n r ) (13) 

p = max(n t , n r ). (14) 


S=[Si s 2 


and 


X = 


Xi 

x 2 


The sum-rate capacity of the MAC in (10) is given by 

C(p) = --supI(X; Y) (15) 

t C 


where the supremum is over all probability distributions on X 
for which 

1) {X I }"2 :1 are independent; 

2) the per-user power constraint 


E 


Tr{XjXj} 


t c Tli P 

< i = l,2,...,n a (16) 

n t 


is satisfied. 


Here, p can be thought of as the total energy per channel use 
available over all users. The particular form of average power 
constraint in (16) allows all users to transmit at the same average 
power per antenna. As reviewed in Section I, C(p) is not known 
in closed form. 

In the next section, we provide nonasymptotic upper and 
lower bounds on the sum-rate capacity C{p) in (15). These 
bounds will be numerically evaluated in Section IV, and con¬ 
clusions be drawn in Section V. 


III. Bounds on Capacity 

Theorem 1 below provides a nonasymptotic upper bound on 
C(p). The upper bound is obtained by dropping the requirement 
that the {Xj}"2. 1 are independent (which enlarges the set of 
distributions over which the maximization in (15) is performed), 
and by bounding the mutual information in (15) using the du¬ 
ality bound [15, Eq. (186)]. The duality approach requires the 
specification of an auxiliary p.d.f. on the channel output Y, 
which, following [8], we choose so that i) Y is isotropic; ii) 
the largest l singular values of Y are distributed as the singular 
values of the noiseless channel output SX, with X following a 
USTM distribution; iii) the remaining n T — l singular values are 
distributed as the singular values of an additive noise matrix of 
appropriate dimension. 

Theorem 1: The sum-rate capacity C{p) in (15) of the 
Rayleigh block-fading MAC (8) is upper-bounded by 

C(p) < u(p) + ^ inf sup g(D, A). (17) 

tc A>0 ])^])n t Xn t 


u{p) = - n r + -p- log ^ 


t c p 

n t 


■log 


P(n r 


n,n r (n r - 
+ + ■ L - ! - 
t c p 

(tc - n t )(n T 


t c 


-l)V(t c -l)P(nt) \ 
V(t c )V(p -1) ) 

- Q(*c - 0 

tc 

— log (p) 


and 


1 

tc 


log (k(p)) 


(18) 


tcp 

+ (t c — n t ) Eg [logdet(G(/ ni + D 2 )+ pl nt )\ 
- n r logdet(J„ t + D 2 ) + X(t c p - Tr{£> 2 }) (19) 


where V(-) is given in (4), G is an n r x n t complex random 
matrix with i.i.d. 1) entries and Jl is the expected value 

of the square of the largest singular value of an n r x (t c — nt) 
complex random matrix with i.i.d. CJ\f(0, 1) entries. Finally, 
with k(p) we denote the probability that the lowest nonzero 
singular value of an n r x n t complex random matrix with i.i.d. 
CAf(Q, tcp/nt) entries is greater than the largest singular value 
of an independent (n r — l) x (t c — l) complex random matrix 
with i.i.d. C7V(0,1) entries. 

Proof: See Appendix B. ■ 

When n t = n r , the sum-rate upper bound (17) can be tight¬ 
ened. This result is given in the following corollary. 

Corollary 2: The sum-rate capacity C{p) in (15) of the 
Rayleigh block-fading MAC (8) for the case n t = n r = n is 
upper-bounded by 

C(p) < u*(p) + inf sup g*(D, A). (20) 

tc X>0 £) gD nxn 


Here, 


and 


u*(p) = -n + — log 
tc 


t c p \ n?_ 
n ) t c p 


■log 


P(n)V(t c ~n) \ 

P(tc) ) 


g*(D, A) = n Tr ^' P ^ — n log det (/„ + D 2 ) 

tcP 

+ A [t c p - Tr{D 2 }) 

(f c - n) ELt det (R k {I n + D 2 )) 
+ IC{D 2 ,tc-n) det(/ rl + D 2 )~ tc+n+1 


( 21 ) 


( 22 ) 


where V(-) is given in (4) and /C(-, •) is given in (6). The matrix 
Rk (A) in (22), which is a function of the nxn diagonal matrix 
A, is defined as follows. Let a\ > a 2 > • • • > a n denote the 
ordered diagonal entries of A. Let the nxn real matrix Pk{A) 
and the (t c — n) x n real matrix Tk, k = 1, 2,..., n be defined 
as follows: 


[Pk(A)]ij 


of fc+1 (log (ai) + ijj{n-k+l)), 


n— 7-f-l 

a,- , 


j = k 


(23) 


1 < * < n, 1 < j < n, j ^ k 














5 


and 


[Tk\ij 


1 

P{n — j + l,t c -n-i)’ 

1 < * < n, 1 < j < n, j ^ k 
if>(t c — i — k + 1) 

/3(n — k + 1, t c — i) ’ 

< l<i<t c — n — 1, j = k (24) 

1, 

i = t c -n , 1 <j<n,j^k 
ip(n — fc + 1), 

i = t c — n, j = k 


where /?(•,•) is defined in the notation section. Finally, let 

Q(A) £ M" x ( tc_n ) and 5 e R( tc_n ) x ( tc_n ) be given by 

[Q(A)]jj = (-ai)- 7-1 "” -4 ', 1 < i < n, 1 < j < t c - n (25) 


and 


Sij = 


(3(i -j + l,t c -n-iy 

(-i y~ j , 


1 < j < i < — n — 1 

l<j<i = t c —n 


0, 


\ < i < j <t c — n. 

(26) 


Then, 


Rk(A) = Pk(A) - Q{A)S~ 1 T k . (27) 

Proof: See Appendix C. ■ 

In Theorem 3 below we provide a lower bound on C (p) . This 
lower bound is obtained by evaluating the mutual information 
in (15) for the MAC-USTM input distribution introduced in 
Section I. Specifically, we set 



with V* uniformly distributed on S(rii,U)- With this choice, 
the power constraint in (16) is satisfied with equality. Similar 
to the bounds developed in [9]—[11], our lower bound relies 
on the Itzykson-Zuber integral [12, Eq. (3.4)]. In order to give 
a compact expression for our lower bound, we shall focus on 
the setup where each user is equipped with a single antenna 
i.e.. Hi = 1, i = 1,2,..., n u . We will address the case when 
the users have multiple antennas at the end of this section. 

Theorem 3: The sum-rate capacity C{p) in (15) of the 
Rayleigh block-fading MAC (8) for the case when all transmit¬ 
ters have a single antenna, i.e., rii = 1, i = 1 , 2 ,..., rz u is 
lower-bounded as follows: 

c( ' >)ilog (^r) + " , ‘ o '’ 

- n r E[log det (/„, + D 2 )] 


+ E 

tog^YY^c-n,))] 


- E 

log ^E d 

det(M(YYt,E)) 

v 

det(/„, + D 2 )™ r /C(E, t c - n t ) 

)\ 


( 29 ) 


Here, D is an n t x n t real diagonal matrix with entries containing 
the singular values of the input matrix X, which is assumed to 
follow the MAC-USTM distribution (28). Furthermore, the n t x 
n t matrix E is defined as E = (D -2 + / n| ) _1 , V(-) is given 
in (4), and is given in (6). Finally, the p x p real matrix 

M(A, B) is defined as follows: 


[M(A,B)]ij 


' 7 (aibj,t c 



p ), 1 < i < n r , 1 < j < nt 

rh + 1 < i < p, 1 < j < n t (30) 
1 < i < n r , n t + 1 < j < p. 


Here, a\ > <22 > • • • > a Ul are the ordered eigenvalues of 
the positive-definite matrix A £ C n ' x n '; similarly bi > b 2 > 

• ■ • > b ril are the ordered eigenvalues of the positive-definite 
matrix B £ C ntXn ‘; the function 7 (-, •) in (30) is given in (7). 

Proof: See Appendix D. ■ 

The lower bound (29) involves expectations that are not 
known in closed form. Hence, we resort to Monte-Carlo meth¬ 
ods for the evaluation of (29). One exception is the two-user 
case n u = 2, for which the expectation over D in (29) admits 
a closed-form integral expression. This result is presented in 
Corollary 4 below. The proof of this corollary exploits that for 
the two-user case, the probability distribution of the eigenvalues 
of XXt can be obtained from the p.d.f. of det(XX^), which 
can be calculated through Bartlett’s decomposition [18, Prop. 
2.1]. This approach does not extend to the case n u > 2or n t > 2. 


Corollary 4: The sum-rate capacity C(p) in (15) of the 
Rayleigh block-fading MAC (8) with n u = 2 and n, = 1 , i = 
1, 2 is lower-bounded as follows: 


C{p) > log 


P(tc - 0 


+ nf c p + log 


Up 

2 


nu) 

j log (^ 2 — a) (1 — a) tc ~ 2 da 
J 0 

E[log(/C(YYt,f c - n r ))] - log(f c - 1) 


- nft c - 1 ) 


- E 


log (j det(M(YY\E(a))) (p 2 -a 2 ) tc ” r 1 da 


(31) 


Here, p = 1 + 2/t c p and 

£(„) = di»g{ 1+ 1 p + a +a) . 1+ 1 M1 ° a) }. (32) 

Furthermore, M(-, •) is given in (30), V(-) is given in (4), and 
/C(-, •) is given in (6). 

Proof: See Appendix E. ■ 

For the case when the users are equipped with multiple an¬ 
tennas, the singular values of the input matrix X are no longer 
distinct, and (29) needs to be further simplified using L’Hopital’s 
rule. The final expression of the resulting lower bound is omitted 
because it is involved. Instead, numerical results are provided in 
Section IV. 


IV. Numerical Results 

Focusing on the setup where the total number of transmit 
antennas n t is equal to the number of receive antennas n r , we 
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Fig. 3. Upper bounds (20) and (33), and lower bound (29) on the sum-rate 
capacity of the MAC (8); n u = 2 single-antenna users, n r = 2, f c € {4, 10}. 



Fig. 4. Upper bounds (20) and (33), and lower bound (29) on the sum-rate 
capacity of the MAC (8); n a = 3 single-antenna users, n r = 3, fc G {6,10}. 


numerically evaluate in this section the upper bound (20) and 
the lower bound (29) on the sum-rate capacity (15). 3 

We also evaluate the i.i.d. Gaussian lower bound developed 
by Rusek et al. [9], which is obtained by setting X = y 7 p/n t G 
in (15), with G having i.i.d. CA/"(0,1) entries. For this choice 
of X, the expectation over D in (29) can be evaluated in closed 
form using the properties of Wishart matrices and the integral 
formula [19, Cor. 2], See [9], [11] for details. 

We also consider the simple capacity upper bound obtained by 
assuming that the receiver has perfect knowledge of the channel 
matrix S in (11). This bound is given by [20] 


C{p) < E 


log det 



(33) 


We start by considering the case of single-antenna users. In 
Figs. 3-5, we depict the upper bounds (20) and (33), and the 
lower bound (29) on the sum-rate capacity of the MAC (8) for 
different values of number of users n u , number of receive anten¬ 
nas n r , and coherence interval t c . For the choice of parameters 
in Figs. 3-5, the upper bound (20) and the lower bound (29) 
characterize the sum-rate capacity accurately already at SNR 
values as low as 10 dB. The tightness of the bounds increases 
as the coherence interval t c or the SNR increases. This last ob¬ 
servation comes as no surprise, since the choice of the auxiliary 
distribution in the derivation of the upper bound (20) is dictated 
by high-SNR considerations. 

The tightness of our bounds implies that the sum-rate capacity 
of the MAC (8) is well-approximated by the capacity of an 
equivalent point-to-point MIMO channel with the same SNR 
and the same total number of transmit and receive antennas. This 
follows because the upper bound is derived under the assump¬ 
tion of perfect cooperation among the users. Furthermore, we 
observe that the perfect-CSI sum-rate capacity upper bound (33) 


^Numerical routines implementing the upper 

bound (20) and lower bound (29) can be downloaded at 
https://github.com/infotheorychalmers/mac_capacity_bounds 



Fig. 5. Upper bounds (20) and (33), and lower bound (29) on the sum-rate 
capacity of the MAC (8); n a = 4 single-antenna users, n r = 4, t c G {8,10}. 


is loose for the relatively small t c values considered in this 
section. 

As shown in Fig. 2 (see Section I) for the case n a = n T = 4 
and t c = 10, the MAC-USTM lower bound is tighter than 
the i.i.d. Gaussian lower bound in the high-SNR regime, al¬ 
though marginally so. The same consideration holds for the case 
n a = n r = 2 and n a = n r = 3 considered in Fig. 3 and Fig. 4. 
Furthermore, the gap between the two bounds gets smaller for t c 
values smaller than 10. 

In Fig. 6, we consider the case when n t = n T = 4andf c = 10, 
but the available 4 antennas are divided unevenly among the 
users. Specifically, we assume n u = 2 (two users) and that the 
first user has one antenna ( n\ = 1) whereas the second user 
has three antennas (712 = 3). Note that the upper bound and the 
i.i.d. Gaussian lower bound depend only on the total number of 
antennas n t and not on the way the antennas are divided among 
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Fig. 6. Upper bound (20), MAC-USTM lower bound (29), and i.i.d. Gaussian 
lower bound [9]; n u = 2 with n\ — 1 and ri 2 = 3; furthermore, n r = 4 and 
t c = 10. 


users. Hence, these two curves coincide with the curves drawn 
in Fig. 1 and Fig. 2. On the contrary, the MAC-USTM bound 
depends on n i and n 2 , and, for the asymmetric antenna scenario 
considered in Fig. 6, is tighter than the i.i.d. Gaussian bound, 
although marginally so. 

V. Conclusion 

We presented finite-SNR upper and lower bounds on the sum- 
rate capacity of Rayleigh block-fading channels for the scenario 
where neither the transmitters nor the receiver have access to 
a priori CSI. The upper bound, which is based on duality, is 
derived under the assumption that the transmitters can cooperate 
perfectly. This transforms the MAC into an equivalent point-to- 
point MIMO channel. We obtain the lower bound by assuming 
that a USTM signal is independently transmitted by each user 
over its available antennas (MAC-USTM). The gap between the 
upper and the lower bounds is less than 8% at 10 dB for the 
case of 4 single-antenna transmitters communicating with a 4- 
antenna receiver over a MAC with coherence interval 10. The 
gap reduces further when the SNR or the coherence interval 
get larger. This implies that the capacity gains obtainable by 
allowing cooperation among transmitters are minimal in this 
case. Our numerical results show also that MAC-USTM yields 
rates that are only marginally larger than the ones obtainable 
using i.i.d. Gaussian inputs. This suggests that i.i.d. Gaussian 
inputs are almost sum-rate capacity optimal for the Rayleigh 
block-fading MACs considered in Section IV. 

Appendix A 
Preliminary Results 

A. An Integral Formula 

The following lemma is useful in the evaluation of integrals 
that involve matrix determinants. We shall need this lemma in 
the proof of Corollary 2. 


Lemma 5: Let o, b £ K., with a < b. Let {/*(•)}?=! and 
i be arbitrary integrable functions over [a, b). Let A = 
A(Xi,X 2 , • ■ ■, x n ) and B = B(x i, £ 2 ,..., x n ) be n x n and 
to x m matrices (m > n) whose entries depend on the scalars 
X \, a:2, • ■ •, x n as follows: 


Aij — 


fi ( x j ) j 


Hu 


Eij — Hi [xj ), 


1 < i < m, 1 < j < n 
1 < i < m, n + 1 < j < m 
1 < i < n, 1 < j < n 


(34) 

(35) 


where c,, are arbitrary scalar real constants. Finally, let E £ 
R mxm be defined as 


Eij — 


f b 

/ fi ( x )9j (x)dx. 1 < i < m, 1 < j < n 
J a 




(36) 

1 < i < m, n + 1 < j < m. 


Then 


J ' J det(A) det(B)dxidx 2 ■ ■ ■ dx n = det(-E'). 

(37) 


a<x n <a; n _i"-<xi <b 


Proof: See [21, Lem. 2]. 


B. Limits of Determinants 

The following lemma, which characterizes the limiting be¬ 
havior of the ratio between the determinant of a certain matrix 
and a Vandermonde determinant, will be needed in the proof of 
Corollary 2. 

Lemma 6: Let A £ C mxm be a positive-definite matrix with 
ordered eigenvalues a\ > > ■ • • > a m . Let C £ C mxm 

be a matrix with entries Cij = fi(a,j), for some differentiable 
functions {/,(•)}£!!i- Then for every integer n < to, and every 
scalar real constant ao, 


lim 

{ a n + l , a n + 2 1 ■ ■ ■ , a m}—> a O 


det(C) 

V(A) 


det (E) 

V(m — n)/C(Ao — ao! n , m — n) 


(38) 


where A 0 = diag{ai, a 2 ,..., a n }, V(-) is given in (4), V(-) is 
given in (5), /C(-, •) is given in (6), and the entries of the m x m 
matrix E are 


Ei. 


fi{af), 1 < i < to, 1 < j < n 

/} m - J )(a 0 ), 1 < i < m, n + 1 < j < m.^ 


Here, f^' Uf) denotes the fcth-order derivative of /) 
Proof: See [22, Lem. 5]. 


C. Expectation of the Log Determinant of a Gaussian 
Quadratic Form 

The following lemma gives a closed-form expression for 
E[logdet(XLX^)] where X has i.i.d. CA/”(0,1) entries and L 
is a certain positive-definite matrix. A closed-form expression 
for the case when the eigenvalues of L are distinct is provided 
in [23, Lem. 2]. Here, we derive a different closed-form expres¬ 
sion, which does not require the eigenvalues to be distinct, and 








appears better suited for numerical evaluations. This lemma is 
used in the proof of Corollary 2. 

Lemma 7: Let X G C nxm , where n < to, be a random 
matrix with i.i.d. CAT(0, 1) entries. Let L G C mxm be a positive- 
definite matrix whose largest n eigenvalues satisfy /i > l 2 > 
••• > In > 1, and whose lowest m — n eigenvalues are equal 
to 1. Then 


E[logdet(XLX t )] 

det(L 0 ) m_n_1 
/C(L 0 - I n ,m - n) 


^det(f4(L 0 )) 

fc =1 


(40) 


Here, L 0 G B nxn is a diagonal matrix, which contains the n 
largest eigenvalues h > l 2 >■■■> l n of L on its diagonal. 
Furthermore, the to x m matrix G(L 0 , E) is defined as follows 

[G(L 0 ,E)\ij = 


exp( 


1 

< 

i <n. 

, 1 <j<n 

m-i 

exp (-ej), 

n 

+ 

1 < i 

< to, 1 < j < n 

(-h) 

j-m 

1 

1 

< 

i < n. 

, n + 1 < j < to 

r(m 
r (i- 

-j + 1 ) 

, n 

+ 

1 <3 

< i < m 

0 , 


n 

+ 

1 < i 

< j < m 


where /C(-, •) is given in (6), L 0 = diag{/i, l 2 , ■ ■ •, l n }, an d 
Rk(-) is defined as in Corollary 2. 

Proof: We first obtain the joint p.d.f. of the eigenvalues of 
S* = XL*Xi, where L* g C mxm is a positive-definite ma¬ 
trix whose largest n eigenvalues coincide with the eigenvalues 
of L , and whose remaining m — n eigenvalues are distinct, i.e., 
Z* + i > (* +2 > ■ • • > 1^. We then derive the p.d.f. of the eigen¬ 
values of S = XLXi by letting {l* +1 , l* l+2 -..., /*,} —> 1. 
Knowledge of this p.d.f. allows us to obtain the moment gen¬ 
erating function of logdet(S) in closed form, from which (40) 
follows by evaluating the derivative of this moment generating 
function at zero. The p.d.f. of S* is [24, Sec. 2] 4 


/s* (S) 

1 det (F(—(L*)-\S)) 

7 , 71 ( 71 — 1)/ 2 v(S')V(-(L*)- 1 ) det(L*) n ’ 


S G C nxn 

(41) 


where F(A, B) is an to x to matrix that depends on the ordered 
eigenvalues {ai}^l 1 and {&i}" =1 of the Hermitian matrices A g 
C mxm and B G C nxn as follows 


[F(A,B)]ii 


„cnbi 


m-j 


1 < i < m, 1 < j < n 
1 < i < m, n + 1 < j < m. 


(42) 


where e\> e 2 > ■ ■ ■ > e n denote the diagonal entries of E. To 
obtain (45), we used Lemma 6, and (46) follows because for all 
A G C nxn , we have that V{—A) = V(A) and that V(A~ l ) = 
V{A)/ det(vl)"- 1 . 

Next, we evaluate the moment generating function of 
log det(S): 

g{t) = E[exp(t log det(S))] = E[det(Es) 4 ] 

= j ••• J det {EYf Es (E)dE 

jjnxn 

det(L 0 ) m " n_1 det (H(t, L 0 )) 

V(n)V(m — n)K.(Lo — I n ,m — n) 

Here, the m x m matrix II (t. Lf) is defined as follows: 

T(t + n — j + 1 )Z* +n_J+1 , 1 <i <n, 1 < j <n 


r(t + to — i + n — j + 1), 

n + l<*<TO, l<j<n 


1 < i < n, n + 1 < j < to 

T(m-j + l) 

r(i - j + 1 ) j ’ 

n+l<j<i<m 

0 , 

n + 1 < i < j < to. 


(51) 


(48) 

(49) 

(50) 


Let Es * denote the n x n real diagonal matrix having the ordered 
eigenvalues of S* on its main diagonal. The p.d.f. of Es* is [26, 
Thm. 3.2] 


To obtain (50), we used the integral formula (37) in Lemma 5. 
To establish (40), we now proceed as follows. Let g'(t) denote 
the first derivative of g(t). Then 


/e s *(£) 

1 det (F(-(L*)~\E))V(E) 
V(n) V(—(L*) -1 ) det(L*) ra 


E G B nxn . (43) 


Finally, let Es denote the n x n real diagonal matrix hav¬ 
ing the ordered eigenvalues of S = XLX^ on its main di¬ 
agonal. We obtain the p.d.f. of Es by computing the limit 

{l* n+1 ,l* n+2 ,...,l* n J-Min (43): 


/e s (£) 




+!>' 


lim 


/E S .(^) 


V{E) _ det (G(L 0 ,E)) 

V(n) det(Lo) n V{m — n)K.(— Lq 1 + / n , m — n) 
det(£ 0 ) m -"- : 1 det(G(L 0 ,E))V(E) Dnx 

V{ri)V{m — n)JC(L 0 — I n , m — n) ’ 


(44) 

(45) 

(46) 


E[logdet(Y)] = lim g'(t) (52) 

_ det(£ 0 ) m -"- 1 ELidet(g fc (L 0 )) 
T(n)V(m — n)IC(L 0 — I n ,m — n) 

where the m x m matrices Hk(L 0 ), k = 1,2,..., n, are given 
by 


[H k (Lo)h = 

f [ff(0,Lo)]y, 

1 < i < m, 1 < j < rn, j k 
r(n — k + l)l™- k+1 (\og{h) + ip(n — k + 1)), 

1 < i < n, j = k 

r(m — i + n — k + \)ip(rn — i + n — k + 1), 

n + 1 < i < m, j = k. 


(54) 


4 There seems to be a typo in the expression given in [24, Sec. 2], The term In (53), we have written the derivative of the determinant of the 

7 , 71 ( 71 — 1)72 should h e i n denominator (see [25, Eq. (50) and Eq. (57)]). matrix H(t, Lf) as a sum of n determinants by using that the 
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matrix H(t , Lq) depends on t only through its first n columns. 
Through algebraic manipulation one can show that 

det (H k (L 0 )) = V(n)V(m - n ) det (R k (L 0 )) (55) 

where Rk{-) is defined in (27). The equality in (55) relates the 
determinant of an m x m matrix to the determinant of a smaller 
n x n matrix. Substituting (55) into (53) we obtain (40). ■ 


Appendix B 
Proof of Theorem 1 

We upper-bound the sum-rate capacity C(p) by enlarging the 
set over which the supremum in (15) is computed. Specifically, 
we drop the assumption that the {Xi}"2. 1 are independent, and 
we substitute (16) with the “global” power constraint 

E[Tr{XX t }] < t c p. (56) 

Let J ip) denote this enlarged set. Then 

C(p) < — supX(X;Y). (57) 

*c J(p) 

Note that the right-hand side (RHS) of (57) is the capacity of an 
n t x n r MIMO Rayleigh bock-fading channel. Hence, it follows 
from [6, Thm. 2] that we can restrict the supremum in (57) to 
input distributions for which X = DQ, where D is an n t x rp 
real diagonal matrix and Q is independent of D, and uniformly 
distributed over the Stiefel manifold <S(n t , f c ). 

Next, we upper-bound the RHS of (57) using the duality 
bound [15, Eq. (186)]. This requires the specification of an 
auxiliary p.d.f. r Y (-) on Y, which we choose as follows. Let 
us denote the singular value decomposition (S VD) of Y as 

Y = UIV. (58) 


Here, U belongs to the set U (n r ) and has real and positive entries 
on the main diagonal, 5 I is an n r x n r diagonal matrix that 
contains the ordered singular values of Y on its main diagonal, 
and V belongs to the Stiefel manifold S(n r ,t c ). Furthermore, 
we choose the auxiliary p.d.f. so that the matrices on the RHS 
of (58) are mutually independent, which leads to 

>v(n = /u(^)/i(S)/v(V)J(E) (59) 


where /u(-), /l(-) and /v(') are the p.d.f. of U, £, and V, 
respectively and 


V(S 2 ) 2 det(S) 2 (^-"d+i 


(60) 


is the Jacobian of the SVD [7, App. A]. Moreover, we assume 
that U and V are uniformly distributed (with respect to the Haar 
measure) over their respective domains, which implies that 


fuAU) 

MV) 


nr(nr-l) 

7T 2 

Hu) 

V(t c — n T )2 nT ir n,tc ” 


(61) 

(62) 


The p.d.f. of £ is chosen as in [8, Sec. V]. Specifically, we 
take the first l singular values of Y to be distributed as the 


singular values of the noiseless channel output SX, with X 
USTM distributed. We take the remaining n r — l singular values 
of Y to be independent of the first l singular values and to 
be distributed as the singular values of an (n r — l) x (f c — l) 
matrix with i.i.d. CAf(0 ,1) entries. Let £i be an l x l diagonal 
matrix containing the first l singular values of Y and let £2 be 
an (n r — l) x (n r — l) matrix containing the remaining singular 
values. Then we have 

r (y y \ /til (^1)/tl 2 (^2) 

1)^2) =- -p-T -• (63) 

K(p) 

Here, 

fc(£ ' ) = piif5r“ I ’ll l l) j 

•det(£i) 2 ( p - z ) +1 V(£?) 2 , Si€D , xI (64) 


where /3 = t c p/n t and 


fz,(S 2 ) 


V(n T - l)V(t c - l) 

•V(£|) 2 , S 2 G O n ’- lxn ’- 1 . (65) 


Furthermore, n(p) is a normalization constant that ensures that 
the singular values are ordered. Specifically, n(p) is the probabil¬ 
ity that the smallest singular value in Li is larger than the largest 
singular values in £ 2 - Using ry(-) in the duality bound [15, Eq. 
(186)] we obtain 

C(p) < j- sup {-E[log(r Y (Y))] - h(Y\X)}. (66) 

Ic J(p) 

Note that the expectation in (66) is not with respect to r Y (•) but 
with respect to the probability distribution on Y induced by the 
input distribution on X through the channel (10). Fix A > 0. 
We can further upper-bound the RHS of (66) by using (56) as 
follows: 

C(p) < ^ inf sup { - E[log(r Y (Y))] - /i(Y|X) 

tc x >°J(p) 

+ A(f c p-E[Tr{XXt}])}. (67) 


Substituting (61), (64), (65), (62), and (60) into (67) and us¬ 
ing [6, Thm. 2] we can rewrite the upper bound (67) as 


C( P )< 


inf sup 
A>0 J(p) 


n T t c + n T n t log(/3) + log(/t(p)) 


+ log ( wrn^D ) 

- n r E[logdet(/„, + D 2 )] + A (t c p 

+ iE[Tr{£ 2 }]+E[Tr{£ 2 }] 

+ (l + t c —p — n r )E[logdet(£ 2 )] 


E 


log 


V(£ 2 ) 


2^2 


V (£ 2 ) 2 V (£ 2 ) 2 


E[Tr{D 2 }]) 


( 68 ) 


We proceed now by upper-bounding some of the terms on the 
RHS of (68). First we bound E [Tr{If}]: 


5 This second requirement ensures that the SVD in (58) is a one-to-one map. E[Tr{Z^}] <jE[Tr{YY^}] = Tl r t c + 7l r E [Tr{D 2 }]. (69) 
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Let a\ > o\ > • ■ • > af M denote the eigenvalues of YY'. To 
bound E[Tr{Z|}], we use the argument given in [7, p. 377], 
which yields 


E[Tr{I^}] = E 


E< 

i=Z+1 


< (n r -/)(fc-0- (70) 


We also have that 


E 


log 


V(£ 2 


V(£?) 2 V( S 2 )2 


=E E E 2 log(< l 2 -^ 2 ) 

i=l j=l+1 

< 2(n r — ()E[logdet(l 2 )]. (71) 


Finally, to bound the term E[logdet(lf)], we start by not¬ 
ing that, given X = X = DQ. the rows of the random 
matrix Y are i.i.d. complex Gaussian random vectors with 
zero mean and covariance matrix I tc + X' X. This means that, 
given X = X = DQ, the matrix Y^Y has the same dis¬ 
tribution as G(/„, + D 2 )G t + HHl where G £ O xn > 
and H £ C n r x(t c -n t ) are i nc j e p enc i en t matrices with i.i.d. 
CA/"(0,1) entries. Let Ai, A 2 ,..., A 1 denote the l eigenvalues of 
G(J„ t + D 2 ')C' (note that this matrix has rank l). Furthermore, 
let fii denote the largest eigenvalue ofHHA, and let Ji = E[/zi]. 
We have 


E [log det (I 2 )] =E 


E lo sK 2 ) 


= E- 


D.Q 


E- 


< Ed 


< Ed 


Y | D, Q 

l 


E, 


G.H 


E lo sK 2 ) 

. 2=1 

l°g(Ai + Ml) 


D,Q 


D 


. 2=1 


E, 


G 


log(Ai + n) 


,i—l 


D 


(72) 

(73) 

(74) 

(75) 


= E[logdet(G(/ n , + D 2 )Gl + Jlln,)] 

- (n r - l)\og(ji). (76) 


Here, (74) follows from Weyl’s theorem [27, Thm. 4.3.1] and 
in (75) we used Jensen’s inequality. Substituting (69), (70), (76) 
into (68) we obtain 


C(p) < u(p) + i inf sup E[p(D, A)] (77) 

^>0 j( p ) 

< u(p) + — inf sup g(D, A) (78) 

t c A>0 £)gI|n|Xnt 

where u(-) and g{-. •) are given in (18) and (19), respectively, 
and where the last step follows by upper-bounding the supre- 
mum over the probability distribution on D with the supremum 
over the set of deterministic diagonal matrices D. 


Appendix C 
Proof of Corollary 2 


When n t = n r = n, the upper bound (68) becomes 

C{p)<j- inf sup |n 2 log(/3) + log ( ^ 
tc A>0 J(p) [_ V V{t c ) 

— nt c + (f c — n) E [log det (YY*) ] 

— n E [log det (/„ + D 2 )] 

+ i E[Tr{YYf}] + \{t c p — E[Tr{D 2 }]) |. (79) 

Differently from the general case treated in Appendix B, in 
the square case some of the terms on the RHS of (79) can be 
computed in closed form. Specifically, 

E[Tr{YY f }] =nf c + nE[Tr{D 2 }]. (80) 

To evaluate E[logdet(YYf)] we use that, given X = X = 
DQ, the rows of Y are i.i.d. complex Gaussian random vectors 
with zero mean and covariance matrix I tc + X f X. This means 
that given X = X = DQ, the matrix YY has the same 
probability distribution as Z(J tc + X t X)Z t , where Z £ C" xtc 
has i.i.d. CAf(0, 1) entries. Hence, we conclude that 


E[logdet(YY t )] 

= E x [E z [log det (Z (4 + Xt X) Z t) ] ] 
det(/ n + D 2 ) 


= E 


\ t c —n—l n 


/C(D 2 , t c — n) 


E det (Rk(I n + D 2 )) 


fe= 1 


(81) 

(82) 


Here, the last step follows from Lemma 7 in Appendix A. Sub¬ 
stituting (80) and (82) into (79) we obtain 

C(p) < u*{p) + E inf sup E[g*(D, A)] (83) 

tc X >°J(p) 


where u*(-) and <?*(-, •) are given in (21) and (22) respectively. 
We conclude the proof by noting that 


E[< 7 *(D, A)] < sup g*(D, A) (84) 

De D nxn 

which yields (20). 


Appendix D 
Proof of Theorem 3 

We lower-bound C(p) by evaluating the mutual information 
on the RHS of (15) for the MAC-USTM distribution (28), which 
yields 

C(p)>Z(X;Y). (85) 

Let X = UDQ denote the singular value decomposition of X, 
where U belongs to the unitary group U(n t ), the diagonal matrix 
D £ B'"' x contains the ordered singular values of X on its 
main diagonal, and Q belongs to the Stiefel manifold S(n t , t c ). 
We next decompose the mutual information as 

T(X; Y) = h(Y) - h( Y|X) (86) 

and evaluate the two terms on the RHS of (86) separately. Since 
Y is conditionally Gaussian given X, the second term in (86) 
can be simplified as 

h(Y|X) = n r E [log det (/„, +D 2 )] +n r t c log(Tre). (87) 
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We next evaluate h(Y). Since X V ~ X for every deterministic 
unitary matrix V G U(t c ). we conclude that Q is uniformly 
distributed over S(n t ,t c ) and independent of D. We also have 
that 


Y = SX + W (88) 

= SUDQ + W (89) 

~ SDQ + W. (90) 

Here, (89) follows from the singular value decomposition of 
X, and (90) holds because S is isotropically distributed, and, 
hence, SU ~ S, which also implies that SU is independent of 
D. The fact that Q is uniformly distributed over S(n t , t c ) and 
independent of D allows us to invoke [11, Prop. 6] and obtain 
a closed form expression for the conditional p.d.f. / Y |d °f Y 
given D. Substituting this closed form expression into 

h( Y) = -E[log(/ Y (Y))] (91) 

= -E[log(E D [/ Y|D (Y|D)])] (92) 

and then substituting (92) and (87) into (86), we obtain (29). 


Appendix E 
Proof of Corollary 4 


We derive in closed form the p.d.f. of the eigenvalues of the 
2x2 matrix XX^, where the two rows Xi G C lxtc and X 2 £ 
C lxtc of the 2 x t c dimensional matrix X are given by 


X,= 



i = 1,2 


(93) 


with Vi € C lxtc and V 2 G C lxtc i.i.d. and uniformly 
distributed on <S(l,t c ) (with respect to the Haar measure). To 
compute the eigenvalues of XX'f, we express it as follows 


XX f 


tcP Vi 
2 V 2 


[Vi 


V 


tcP [ 1 VrV^' 

2 [v 2 v] 1 


(94) 

(95) 


Let a denote the absolute value of the scalar V^Vj. It follows 
from (95) that the eigenvalues of XX^ are t c p( 1 ±a)/2. Hence, 
in order to determine the p.d.f. of the eigenvalues of XX^, it is 
sufficient to obtain the p.d.f. of a. Since 

det(XXt) = (1 -a 2 ) (96) 

Bartlett’s decomposition (see [18, Prop. 2.1]) implies that a 2 ~ 
Beta(l,£ c — 1)- After some mathematical manipulations, we 
obtain that 


f a (x) = (t c - l)2a;(l - x 2 ) tc 2 , 0 < x < 1. (97) 

Using (97) in (29) together with the expression for the eigenval¬ 
ues of XXt, we obtain (31). 
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